Linear magnetization dependence of the intrinsic anomalous Hall effect 
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The anomalous Hall effect is investigated experimentally and theoretically for ferromagnetic thin 
films of MnsGe3. We have separated the intrinsic and extrinsic contributions to the experimental 
anomalous Hall effect, and calculated the intrinsic anomalous Hall conductivity from the Berry cur- 
vature of the Bloch states using first-principles methods. The intrinsic anomalous Hall conductivity 
depends linearly on the magnetization, which can be understood from the long wavelength fluctu- 
ations of the spin orientation at finite temperatures. The quantitative agreement between theory 
and experiment is remarkably good, not only near K, but also at finite temperatures, up to about 
~240 K (0.8 T c ) 

PACS numbers: 75.47.-m, 71.15.-m 



The appearance of a transverse voltage or electric field 
E y in a metal or semiconductor in response to a longitu- 
dinal electric current j x and perpendicular magnetic field 
B is known as the Hall effect 0, 0, . In nonmagnetic 
materials, this transverse voltage arises from a deflection 
of charge carriers by the Lorentz force j x xB, resulting 
in a Hall resistivity p = E y /j x that is proportional to 
the magnetic field for weak fields. Only one year after 
his discovery of this "ordinary Hall effect", Edwin Hall 
found that the Hall resistivity in ferromagnetic metals 
acquires an extra term which depends on the magneti- 
zation of the samples 0]- Subsequent studies Q found 
that this "anomalous" term is proportional to the spon- 
taneous magnetization M. Empirically, one finds 

Ph = Poh + Pah = RoB + RA-kM (1) 

where pon is the ordinary Hall resistivity due to the 
Lorentz force in a perpendicular magnetic field B, pah 
the anomalous Hall resistivity, Rq the ordinary Hall co- 
efficient, and R s the anomalous or spontaneous Hall co- 
efficient. 

Early theoretical interpretations of the anomalous Hall 
effect (AHE) pointed toward asymmetric scattering of 
the spin polarized charge carriers in the presence of spin- 
orbit coupling. Such scattering mechanisms can explain 
most of the qualitative features of the AHE observed in 
experiments, including the linear Q or quadratic [j| cor- 
relation with the longitudinal resistivity p XXl i.e., 

Pah = ap xx + bp 2 xx (2) 

The two terms on the right are traditionally known as 
the skew scattering Q and side jump [j| contributions, 
respectively. However, quantitative agreement between 
scattering theories and experiment remained largely un- 
settled, in part because the scattering potentials are un- 
known. In recent years, inspired by the new insight on 
Berry phase effects on Bloch electrons a number of 



groups have evaluated the intrinsic anomalous Hall con- 
ductivity (AHC) for ferromagnetic semiconductors 0, 
transition metals || and oxides 0, using first princi- 
ples calculations. The intrinsic contribution can be quite 
large and for the first time, quantitative agreement be- 
tween theory and experiment appears quite reasonable. 
Interestingly, this intrinsic effect does not take into ac- 
count scattering by impurities or phonons. It stems from 
spin-orbit coupling in the crystal band structure, a mech- 
anism originally due to Karplus and Luttinger [lfj. 

In this Letter, we examine the magnetization depen- 
dence of the AHE in light of the new theory. We have 
performed measurements on single-crystal films of fer- 
romagnetic Mn5Ge3, extracted the intrinsic AHC, and 
found that this intrinsic conductivity is linear in the mag- 
netization over a wide range of temperatures. This is 
puzzling. Although a linear dependence on magnetiza- 
tion comes naturally from the scattering mechanisms, 
the intrinsic mechanism often yields strongly nonlinear 
behavior 0, 0, 0- We solve this puzzle by invoking 
the spin-fluctuation picture for finite temperature mag- 
netism. First principles calculations for the AHC based 
on this picture yield almost perfect comparison with ex- 
periment. 

Mn5Ge3(0001) thin films were grown on Ge(lll) by 
solid-phase epitaxy, following the procedures in Ref. 11, 
and by co-deposition of Mn and Ge on Ge buffered 
GaAs(lll). Both types of film exhibit very similar trans- 
port properties at low temperature but a GaAs substrate 
is desirable for transport measurements above ^200 K 
because it minimizes the parallel conductance through 
the substrate. For the co-deposition experiment, we 
adopted the following procedure. A GaAs(lll) substrate 
was annealed in ultrahigh vacuum at about 600°C. A 35- 
nm-thick Ge buffer layer was subsequently grown on top 
of GaAs(lll) at 400°C. Next, we co-deposited Mn and 
Ge in a 5:3 atomic ratio onto the Ge(lll)-c(2x8) buffer 
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FIG. 1: (Color online) (a) Temperature-dependent sponta- 
neous magnetization M. The upper inset shows the RHEED 
pattern of a MnsGe3 thin film grown on Ge/GaAs(lll). The 
lower inset shows the temperature dependence of the inverse 
magnetic susceptibility, measured at 0.7 T. The green line is 
a linear fit and extrapolates to Tc = 298 K. (b) M versus T 2 . 
The green line is a linear fit, showing that AM(T) /M(0) <x T 2 
up to ~220 K. 

layer substrate at 330° C. The growth rate was ~6A/ mm, 
and the typical film thickness is 45 nm. Sharp reflection 
high energy electron diffraction (RHEED) patterns con- 
firm the epitaxial growth reported in Ref. 11. Magnetic 
properties were measured with a superconducting quan- 
tum interference device magnetometer. The longitudinal 
resistivity was measured using the standard four-point- 
probe technique. Hall measurements were performed fol- 
lowing the procedures descibed in Ref. 12. 

Fig. 1 shows the spontaneous magnetization, obtained 
by extrapolating the high field part of the magnetization 
isotherms M(H) to zero internal field, after correcting 
for the linear diamagnetic background signal. The satu- 
ration magnetization at 2.5 K is 2.6 ± 0.2 pg per Mn, in 
excellent agreement with our previous results from solid- 
phase epitaxy The susceptibility x — M/B follows 
the Curie- Weiss law (Fig. 1(a): inset), and the Curie 
temperature Tc determined from the linear fit is 298 ± 
3 K. Our temperature range (2 - 400 K) is not low enough 
to resolve the T 3 / 2 dependence in the spontaneous mag- 
netization expected from independent spin-wave excita- 
tions. Instead, we find a near perfect T 2 dependence up 
to about 220 K, as shown in Fig. 1(b). The T 2 fall of the 
magnetization arises from long wavelength, low frequency 
fluctuations of the magnetization, or multiple excitations 
of the interacting spin waves [l3j . 

Fig. 2 shows the Hall resistivity at various temper- 
atures. The anomalous Hall component is well charac- 
terized by the "knee" profile below Tc and can still be 
identified at temperatures as low as 2.5 K (Fig. 2a). Rq 
can in principle be extracted from the high field slope of 
the Hall isotherm [14|. The slope changes from negative 
to positive at about 180 K (Fig. 2(b)), indicating a sign 
change of i?o Pah is obtained by extrapolating the 
high field slope to zero internal field, and is shown in Fig. 
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FIG. 2: (Color online) Field-dependent Hall resistivity pu at 
various temperatures 

butions. In light of the scmiclassical transport theory, 
the electrical current can be written as |8j 

(3) 

where Sl n is the Berry curvature of the Bloch state de- 
fined by ft„(k) = -Im(V k u nk |x| V k w nk ). The func- 
tion u nk is the part of the Bloch wave function (with 
band index n) that is periodic in the lattice. f n is the 
equilibrium Fermi-Dirac distribution function, and Sf n 
is a shift proportional to the electric field and relaxation 
time. The first term is an intrinsic anomalous current 
and originates from the Berry curvature correction to the 
group velocity of a Bloch electron. This intrinsic contri- 
bution is independent of scattering and should lead to a 
quadratic dependence, pah oc p xx . The second term nor- 
mally represents the longitudinal current but it will have 
a transverse component or Hall current in the presence 
of skew scattering. Because Sf n is proportional to the 
relaxation time, the skew scattering contribution is pro- 
portional to p xx . To separate the intrinsic and extrinsic 
contributions to the AHE, we thus write 

Pah = a(M)p xx + b(M)p 2 xx (4) 

where the functions a(M) and b(M) generally depend 
on the magnetization. In principle, it is impossible 
to uniquely separate the intrinsic and extrinsic contri- 
butions from temperature or field-dependent measure- 
ments on a single sample because temperature changes 
M, a{M), b(M), and p xx simultaneously. However, the 
skew scattering contribution a(M) is usually linear in 
magnetization 15]. Accordingly, a(M) can be obtained 
by plotting pAn/Mp xx versus p xx , as shown in Fig. 3(b). 
This plot is linear below p xx « 110 pQcm (or, equiv- 
alently T sa 220 K). The intrinsic contribution to the 
AHE is obtained by subtracting a(M)p xx from the ex- 
perimental pah- The intrinsic anomalous Hall resistivity 
is also indicated in Fig. 3(a). The AHC is expressed as 
&ah = Pah I Pxx E3j so we immediately identify b(M) 
as the intrinsic AHC, uiah- The constant slope in Fig. 
3(b) indicates that b(M)/M is constant, hence the ctjah 
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FIG. 3: (Color online) (a) pah as a function of temperature 
before (black line) and after (red line) subtracting the skew 
scattering contribution. p xx is also shown, (b) pah / {M p xx ) 
as a function of p xx before (black line) and after (red line) 
subtracting the skew scattering term. The green line is a 
guide to the eyes crossing the origin. 



is proportional to M, i.e. 



aiAH oc M 



(5) 



This relation is the key result of our paper. Experimen- 
tally, this relation is valid between 1.7 and 2.6 [Ib as 
shown by the linear ht in Fig. 4(c). It spans an amaz- 
ingly broad temperature interval of about 240 K or 0.8 
Tc- ctiah extrapolates to 860 Sl^cm" 1 at K. Note 
that the uncorrected a ah shows a non-monotonic de- 
pendence on M. 

Based on the semi-classical transport theory |(J, the 
intrinsic anomalous Hall conductivity (AHC) can be ex- 
pressed as a sum of the Berry curvatures Q : 



&IAH = — -r- 

n 



BZ 



d 3 k 



n z (k) 



(6) 



where ft z (k) = J2 n /n^n( k )- The intrinsic AHC can be 
evaluated from the electronic band structure. Following 
the procedures in Ref. 8, we first obtained a fully con- 
verged ground state by using the full-potential linearized 
augmented plane-wave method |l6| with generalized gra- 
dient approximation fl7j | for the exchange-correlation po- 
tential. In this calculation, we sampled 1,000 fc-points in 
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FIG. 4: (Color online) (a) Calculated oiah versus the effec- 
tive spin-orbit coupling strength £ relative to the real value 
£o for MnsGe3. (b) Calculated ctiah versus the magnitude 
of magnetization for different relaxation time r. (c) Experi- 
mental AHC as a function of M before (black dots) and after 
(red dots) subtracting the skew scattering (S.S.) contribution. 
The red line is linear fit. Calculated oiah values versus z- 
component of the magnetization M z are connected by the blue 
line. 



potentials we obtain a fully converged (Jiah by summing 
the Berry curvature using a much larger set of fc-points. 
The final converged ujah value is 964 il^cm -1 at T = 
K, in excellent agreement with the experimental extrap- 
olation of 860 £l~ 1 cm~ 1 . 

Next, we explore various mechanisms that could lead 
to a linear magnetization dependence of uiah for T > 
K. A linear dependence on the magnetization is often ob- 
served and is automatically implied by treating spin orbit 
coupling as a linear perturbation in scattering theory |l5| . 
On the other hand, uiah of e.g. Mn doped GaAs 7] and 
iron is nonlinear in spin-orbit coupling, and it appears 
difficult to obtain (J\ah oc M. More dramatically, in the 
case of SrRuOs Fang et al. also found non-monotonic 
or even spiky dependences on the magnetization, which 
they attributed to "magnetic monopoles" in momentum 
space 0- To check whether oiah is linear in the spin- 
orbit coupling strength, we artificially changed the speed 
of light c. This is equivalent to changing the spin-orbit 
coupling strength £ because £ oc c~ 2 J8j. Fig. 4(a) shows 
that a i ah increases as £ increases, but it is clearly non- 
linear in £. This nonlinearity means that spin-orbit in- 
teraction cannot be treated perturbatively for MnsGea. 
It also implies that <jiah oc M cannot be obtained from 
a perturbative analysis of spin-orbit coupling. 
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Next, we investigate whether o~iah oc M can be ob- 
tained by varying the spin-polarization or exchange split- 
ting of the spin-up and spin-down bands of the Stoner 
model, as was done in Ref. 9. After obtaining self- 
consistent electronic charge density, we first tuned the 
exchange potential, then recalculated the wave functions 
and band structure (here no self-consistent calculations 
is needed and desired), and finally obtained ctiah- The 
result is shown in Fig. 4(b) for various relaxation times 
t. The qualitative and quantitative disagreement with 
experiment is evident. 

We now surmise that the linear dependence on the 
magnetization for o-jah, can be quantitatively accounted 
for by the long-wavelength, low-frequency fluctuations of 
the spin orientation at finite temperatures. The exper- 
iment shows that o~iah is linear in M as long as the 
magnetization falls off quadratically (2 K < T < 240 K) . 
In this temperature interval, the magnitude of magneti- 
zation stays constant but the local magnetization rotates 
away from the z-axis. We assume that the typical wave- 
length is much larger than the mean free path, so that we 
may employ the local approximation in the calculation of 
the conductivity tensor. In other words, we assume that 
the conductivity tensor is defined in each region and that 
the local conductivity can be calculated by considering 
an artificial infinite system with a uniform magnetiza- 
tion given by the local magnetization orientation in that 
region. 

Define the local spin orientation by the polar angles 
(6, ip). The local intrinsic anomalous Hall conductivity 
is calculated as follows: (i) First, the self-consistent den- 
sities for both spin-up and spin-down electrons are ob- 
tained, (ii) Next, the spin quantization axis is rotated 
away from the z-axis to a direction specified by (#,<£>). 
(iii) The Kohn-Sham orbitals are then recalculated in the 
presence of spin-orbit coupling, which arc in the form of 
Bloch spinors. (iv) Finally, the Berry curvatures of the 
occupied Bloch spinors are calculated and summed over 
to obtain the intrinsic Hall conductivity, as was done in 
Ref. [8]. 

The observed intrinsic anomalous Hall conductivity 
should be a spatial average of the local values. The result 
is averaged over the azimuth angle ip and found to be pro- 
portional to the z-component of the local magnetization 
to a very good approximation. This is shown in Fig. 4(c) 
by the open circles. Because of this near-perfect linearity, 
we expect that further averaging over the small spread of 
the local z-components will not change our result. The 
slope as well as numerical magnitude of the theoretical 
data are in almost perfect agreement with experiment, 
provided that the extrinsic skew scattering term is prop- 
erly subtracted. The close agreement also suggests that a 
possible side jump contribution to the AHC (which also 
has pah oc Pxx) must be very small or negligible. Theory 
and experiment deviate in the vicinity of Tq where the 
interacting spin waves produce a different fall off of the 



magnetization. 

In conclusion, the AHC of ferromagnetic MnsGe3 thin 
films has a large intrinsic contribution that varies lin- 
early with the magnetization. The agreement between 
low temperature measurements and T = K ab initio 
calculations is already compelling in its own right. More- 
over, the temperature dependence of the intrinsic AHC 
can be quantitatively accounted for by long wavelength 
fluctuations of the spin orientation up to 0.8 Tc or 240 
K. 
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